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ABSTRACT 


The  specialization  of  the  simplex  algorithm  for  the  solution  of 
generalized  network  flow  problems  rests  on  the  fact  that  a  basis  for 
the  problem  may  be  represented  graphically  as  a  spanning  forest  in 
which  each  component  is  either  a  one-tree  or  a  rooted  tree.  The 
design  of  a  specialized  algorithm  for  efficient  solution  of  generalized 
network  problems  necessarily  depends  on  data  structures  chosen  to 
represent  the  basis.  This  paper  presents  the  design  and  detailed 
algorithmic  specification  of  the  primal  simplex  algorithm  for  such 
problems.  Computational  testing  to  determine  the  overhead  required 
by  generalized  network  data  structures  over  pure  network  data  structures 
indicates  that  generalized  network  algorithms  are  on  the  order  of  2.5 
to  3.5  times  slower  than  pure  network  algorithms.  Computational  testing 
with  generalized  network  problems  with  up  to  1000  nodes  and  7000  arcs 

establishes  the  suitability  of  the  data-structures  for  efficient 

;  ^  ’ 

implementation  of  primal  simplex  calculations.  ^  ' 
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I.  Introduction 

The  generalized  network  probl^  is  a  specialized  linear  progranning 
problem  which  may  be  defined  as 

min  cx 

s.t.  Ax  =  r 

0  ^  X  <  u 

where  A  is  a  nxm  node-arc  incidence  matrix  with  multipliers  on  the  arcs. 

That  is,  each  column,  aj^,  has  at  most  two  non-zero  entries,  in  row  i  and 
g|^  in  row  j.  The  problem  may  be  easily  scaled  so  that  one  of  the  entries, 
say  dj^  is  a  +  1.  If  gj^  is  zero,  then  the  variable  corresponds  to  a  loop. 
Otherwise,  the  column  aj^  represents  an  arc  from  node  i  to  node  j  where  dj^ 
units  leave  node  i,  producing  g^^  units  at  node  j.  Without  loss  of  generality, 
we  assume  that  the  matrix  A  corresponds  to  a  connected  network ,  G(N,E),  where 
the  rows  correspond  to  the  node  set  H  of  the  network,  and  the  variables 
correspond  to  the  arc  set  E  of  the  network.  Further,  we  assume  that  the 
matrix  A  is  of  full  row  rank.  It  can  be  shown  that  if  the  rank  is  n-1,  then 
the  problem  corresponds  to  a  pure  network  problem  [4]. 

The  specialization  of  the  simplex  algorithm  for  the  solution  cf  such 
problems  rests  on  the  fact  that  the  basis  for  a  generalized  network  problem 
corresponds  to  a  graphical  structure  referred  to  as  a  spanning  forest.  Much 
in  the  same  manner  that  the  development  of  specialized  solution  software  for 
the  pure  network  problem  rests  on  data  structures  for  the  representation  of  a 
spanning  tree,  specialized  solution  software  for  solution  of  the  generalized 
network  problem  rests  on  the  representation  of  a  spanning  forest. 

The  design  of  a  specialization  of  the  simplex  method  for  the  solution 
of  generalized  network  problems  is  necessarily  dependent  on  the  data  structure 


chosen  to  represent  the  basis.  The  basis  for  a  generalized  network  flow 
problem  has  the  following  form: 


where  each  of  the  t  submatrices  8^.  has  the  form: 


(a)  1-tree  with  self  loop  '  (b)  1-tree  with  loop 


The  basis  corresponds  to  a  spanning  forest  on  the  node  set  N  in  which 
each  matrix  corresponds  to  a  l-tree. 

A  1-tree  on  a  node  set  N^.  with  cardinality  n^.  consists  of  a  spanning 
tree  on  the  node  set  together  with  one  additional  arc,  thus  forming  exactly 
one  cycle.  The  cycle  may  correspond  to  a  self  loop  as  in  (a)-  This  is  the 
main  distinction  between  the  generalized  network  flow  problem  and  the  pure 
network  problem.  Thus  a  data  structure  which  represents  a  spanning  forest 
should  be  able  to  distinguish  between  subsets  of  nodes.  That  is,  distinguish 
between  the  various  1-trees  in  the  basis.  Further,  the  data  structure  should 
facilitate  the  three  essential  operations  of  a  primal  simplex  pivot: 

(1)  pricing,  or  the  determination  of  an  entering  variable,  C2)  the  ratio  test 


operation  which  finds  the  representation  of  the  incoming  variable  with  respect 
to  the  current  basis  and  then  finds  the  leaving  variable  and  (3)  the  update  of 
the  dua*!  variables  to  effect  the  pricing  operation. 

II .  Data  Structure 

An  examination  of  the  requirements  for  effecting  the  primal  simplex 
operations  makes  it  clear  that  the  pricing  operation  is  trivial  once  given 
the  dual  variables  at  any  iteration.  The  column  uodate  or  the  ratio  test 
operations  are  more  involved.  The  entering  arc  may  connect  two  nodes  in  the 
same  1-tree,  or  it  may  have  one  node  in  one  1-tree  and  the  other  in  another 
1-tree.  There  are  several  possibilities  for  the  leaving  arc: 

(a)  The  entering  arc  is  in  a  1-tree,  and  the  leaving  arc  is  in  the  loop  of 
the  1-tree.  Therefore,  a  new  loop  will  be  formed,  and  the  previous 
loop  will  be  broken. 

(b)  The  entering  arc  is  in  a  1-tree,  and  the  leaving  arc  is  tn  the  loop 
formed  by  the  entering  arc.  Therefore,  there  is  no  new  loop  formed 
and  no  old  loop  broken.  This  is  similar  to  a  pure  network  pivot. 

(c)  The  entering  arc  is  in  a  1-tree.  The  leaving  arc  is  neither  in  the 
old  loop  nor  in  the  new  loop  formed  by  the  entering  arc.  Therefore, 
the  old  loop  remains  intact,  and  a  new  loop  is  formed.  Thus,  the  old 
1-tree  is  split  into  two  separate  1- trees. 

(d)  The  entering  arc  connects  two  1-trees,  and  the  leaving  arc  is  a  loop 
arc.  Thus,  the  two  1  trees  are  coalesced  into  a  single  1-tree. 

Ce)  The  entering  arc  connects  two  1- trees,  and  the  leaving  arc  is  not  a 
loop  arc.  Thus,  after  the  pivot,  there  are  still  two  1-trees  with 
the  old  loops  intact. 

In  order  to  effect  pivots  of  the  type  summarized  in  (a)  -  (e),  it  is 
important  to  distinguish  between  loop  arcs  of  a  1-tree  and  to  be  able  to 


distinguish  between  the  various  1-trees  of  a  forest.  In  general,  since  the 
number  of  1-trees  in  a  basis  can  vary,  a  data  structure  should  be  able  to 
effect  changes  in  spanning  forests  efficiently. 

The  data  structure  used  for  storing  the  spanning  forest  corresponding 
to  the  basis  of  a  generalized  network  problem  consists  of  the  following 
node-length  arrays; 

p(i)  The  predecessor  of  a  node  i, 

b(i)  The  brother  of  a  node  i,  (0  for  looo  nodes), 

s(i)  The  successor  of  a  node  i, 

l(i)  The  level  of  a  node  i.  For  a  node  which  is  not  part  of  a  loop, 
this  gives  the  distance  from  the  nearest  node  of  the  loop.  For 
a  node  which  is  part  of  the  loop,  this  is  either  0  or  a  negative 
number.  Zero  for  a  self-loop.  Negative  for  a  loop  node  of  a 
1-tree  where  the  absolute  value  gives  the  number  of  the  particular 
1-tree,  i.e.  the  name  given  to  the  1-tree. 
a(i)  The  arc  connecting  node  i  to  its  predecessor.  This  number  is 
positive  if  the  arc  is  from  the  node  to  its  predecessor  and 
negative  otherwise. 

In  Figure  1,  this  data  structure  is  illustrated  for  a  spanning  forest 
on  16  nodes  with  two  1-trees.  One  of  the  1- trees  contains  a  self-loop  while 
the  other  contains  a  loop  with  three  arcs.  The  brother  and  successor  labels 
allow  a  traversal  of  the  1-trees.  While  a  preorder  traversal  [  ]  has  been 
shown  to  be  suitable  for  a  tree,  it  cannot  determine  the  distinction  between 
various  subtrees  of  a  1-tree  rooted  on  any  one  of  the  loop-nodes.  A  traversal 
of  a  1-tree  is  possible  using  the  brother  and  successor  labels. 

The  spanning  1-tree  can  be  regarded  as  several  rooted  trees  together  in  a 
loop.  If  this  loop  is  ignored,  then  every  non-loop  node  has  a  unique  path  to  the 
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nearest  loop  node  which  is  its  root.  The  predecessor  of  a  non-loop  node  is 
the  next  node  on  the  path  to  this  root.  For  a  loop  node,  the  predecessor  is 
the  next  node  on  the  loop  so  that  a  trace  of  the  loop  can  be  effected  by 
using  the  predecessor.  For  a  self-loop,  the  predecessor  of  the  root  node  is 
itself.  Thus,  the  predecessor  allows  for  a  traversal  from  any  node  to  the 
root  and  the  entire  loop.  For  examole,  in  Figure  1,  by  making  use  of  the 
predecessor,  the  traversal  from  node  5  down  to  node  1  can  be  effected  and 
further,  the  traversal  around  the  loop  from  node  1  to  node  4  to  node  7  and 
back  to  node  1. 

Apart  from  the  above  mentioned  arrays,  the  following  node  length  arrays 
are  used  in  addition; 

f(i)  the  flow  on  arc  a(i), 

v(i)  the  dual  variable  for  node  i ,. 

y(i)  the  updated  column  representation  coefficient. 

Thus,  a  total  of  8-node  length  arrays  are  used  together  with  at  most 
half-node  length  arrays  for  the  loop  factors.  There  are  only  four  arc  length 
arrays  necessary  for  storing  the  arc  data.  These  are  the  to  node  node(x), 
the  multiplier  or  gain  g(x),  the  cost  c(x)  and  the  capacity  cap(x). 

III.  Ratio  Test 

The  ratio  test  operation  consists  of  evaluating  the  vector  y  such 

that 

By  = 

Since  a.  *  d.  e.  -  e.  ,  where  e.  is  the  i^"  unit  vector,  and  because  of  the 
block  diagonal  structure  of  the  basis,  B,  the  non-zero  entries  of  the  vector 
aj^  may  occur  in  the  same  block  of  8  or  in  two  different  blocks  of  B.  In  the 
primal  simplex  specialization  for  the  pure  network  problem,  the  ratio  test 
can  be  performed  as  the  vector  y  is  evaluated  since  the  entries  in  the 


y  vector  are  either  +1,  -1,  or  0.  Thus,  only  a  single  pass  of  the  relevant 
arcs  in  the  basis  which  produce  non-zero  entries  in  the  y  vector  is  required 
For  the  specialization  of  the  primal  simplex  algorithm  for  the  generalized 
network  problem,  it  is  desirable  to  evaluate  the  y  vector  in  a  manner  such 
that  the  analogous  number  of  passes  is  a  minimum.  Two  cases  arise  immediately 

(1)  The  non-zero  entries  of  the  entering  arc  occur  in  two  different 
blocks,  i.e.,  two  1-trees. 

(2)  The  non-zero  entries  of  the  entering  arc  occur  in  the  same  block 
of  the  basis,  i.e.,  one  1-tree. 

By  viewing  the  operation  as 
By  *  By"  =  d|^e,  -  gi^ej  , 

the  essential  operation  to  be  performed  is  the  solution  of  the  following 
system  of  equations,  the  system  of  equations  for  a  single  block  of  the 
basis  B.  A  single  block  corresponds  to  a  spanning  1-tree.  • 

For  simplicity,  we  consider  only  the  rows  and  columns  of  the  block 
which  correspond  to  non-zero  entries  in  the  y  vector  to  be  evaluated. 

The  system  of  equations  then  takes  the  form: 


Note  that  one  of  the  entries  in  each  column,  either  d^  or  g^  will  be  a  +1, 
and  all  the  entries  d^  are  non-zero. 
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This  corresponds  to. the  graph  shown  in  the  figure  2.  The  incidence 
node  for  the  entering  arc  corresponds  to  the  node  for  row  1  or  the  above 
system.  We  shall  assume  that  the  nodes  are  numbered  in  the  order  of  the 
rows.  There  are  p-1  nodes  until  the  first  loop  node,  node  p,  is  reached, 
and  nodes  p  through  q  are  all  in  the  loop.  Note  that  if  9q  is  zero,  then 
the  node  q  has  a  self-1 ooo,  and  o  =  q.  The  solution  to  the  above  system  is 
obtained  as: 

Yl  =  1/dj 

Yj  ■=  i  =  2,3,...,p-l,p+l . q 

'  ^-1  •  '  S)] 

wnere  F  =  (-9„)  •  •  •  (-9-)  /  (d  . ..d  )  and  G  =  1  -  F,  and  for  self-loop  F  =  0, 

F  is  called  the  loop  factor.  This  calculation  can  be  performed  on  the  1-tree 
by  making  use  of  the  suggested  data  structure.  Tracing  the  path  down  to  the 
loop  is  effected  by  making  use  of  the  predecessor  index.  The  first  node  on 
the  loop  is  identified  by  making  use  of  the  level  index.  Note  that  the  loop 
factor  F,  as  well  as  G,  can  be  obtained  when  the  1-tree  is  constructed.  Thus, 
it  is  not  necessary  to  trace  the  loop  twice  to  effect  the  above  solution  of 
equations.  The  y  vector  corresponding  to  any  entering  arc  can  be  obtained  by 
evaluating  y'  and  y"  from 

By'  • 

By*  - 


and  then  setting,  y  *  dj^'  9^^“* 


Depending  on  whether  the  entering  arc  is  incident  to  two  blocks  or  one 
block  the  ratio  test  and  the  y  vector  evaluation  can  be  performed  simultaneously 
or  simultaneously  in  part  only.  If  there  are  two  blocks  then,  the  ratio  test 


can  be  performed  as  the  vectors  y'  and  y“  are  evaluated.  This  is  because 
no  basic  arc  can  appear  in  both  y*  and  y“.  Thus  no  superfluous  calculation 
is  performed. 

If  there  is  only  one  block  then,  arcs  on  the  path  from  node  i  to  the 
loop  and  arcs  on  the  path  from  node  j  to  the  loop  may  in  fact  overlap.  It  is 
desirable  to  arrange  the  calculations  so  that  the  ratio  test  can  be  performed 
as  the  y  vector  is  evaluated.  In  order  to  perform  this  operation  efficiently 
we  define  a  node  which  is  called  the  junction  of  the  paths  to  the  root.  The 
junction  is  defined  to  be  the  common  node  with  highest  level  on  the  two  paths. 
It  may  not  always  exist.  The  following  algorithm  summarizes  the  ratio  test 
operation. 

Algorithm  1.  Ratio  Test 

0.  [Initialize] 

a.  [set  y]  y(x.).  :*  0  x  =  1,  2,  ...,  m 

0 

b.  [set  junction  flag]  junction  :»  0 

c.  [increase  or  decrease  flow  on  arc  k]  5  :»  sign(cap(k)) 

d.  [set  leaving  arc]  u  :*  0 

e.  [set  A  to  bound]  A  H  cap(k)  II 

1.  Travel  from  node  1  to  its  root  node,  where  the  root  node  is  the  loop 
node  that  has  the  nearest  distance  from  the  node. 

a.  [initialize]  T  :»  1  x  ;»  1 

b.  [root  node  ?]  if  T{x)  s  0  ,  then  T1  :•  T 

[store  the  root  node  of  1]  R1  :*  x 

go  to  2 


c.  [evaluate  y(x)] 


if  a(x)  >  0  ,  then  y(x)  :=  T  +  y(x) 

T  :=  -T  •  g[a(x)] 
go  to  d. 

otherwise  T  :=  T/g[.a(x)] 

y(x)  :=  T  +  y(x) 

T  :=  -T 
go  to  d. 

d.  Determine  maximum  allowable  flow  change 

d,l.  [increase  flow  ?]  if  5  •  y(x)  >  0  then  go  to  d.3 

d.2.  [increase  flow]  if  capT  II  a(x)  il  ]  -  f(x)  >  ^ 

-6  •  y(x) 

then  go  to  d.4. 

otherwise  A  :=  cap[  II  (a(x))  II  ]  ••  f(x) 

-6  •  y(x) 

u  :=  X 

go  to  d.4. 

d.3.  [decrease  flow]  if  f(x)  >  A 

5  •  y(x) 

then  go  to  d.4. 

otherwise  A  :=  f(x)/(5  •  y(x)) 

u  X 

go  to  d.4. 


d.4.  [move  to  next  node]  x  :■  p(x) 
d.5.  go  to  b. 
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Algorithm  1.  continued. 

2.  Travel  from  node  j  to  its  root  node 

a.  [initialize],  T  :=  g(k) 

X  :=  j 

u  1  :  =  u  ,  u  :  =  0 

b.  [root  node  ?]  if  l(x)  <  0  then  T2  :=  T 

R2  :=  X 
go  to  3. 

c.  [junction  node?]  if  junction  =  0  and  y(x)  f  0 

then  junction  =  x  , 
go  to  e. 

d.  [evaluate  y(x)  ,  determine  maximum  allowable  flow  change] 

the  same  as  l.c  and  l.d  ,  except,  change  L  :=  1  to  L  2 

u  u 

e.  [redo  ratio  test  along  the  path  from  node  i  to  node  junction  ?] 
if  >  0  and  1  (junction)  >  1  (u^) 

then  go  to  f. 
otherwise  go  to  d. 

f.  [redo  ratio  test  from  i  to  junction] 
f.l.  [initialize]  x  :=  i 

f.2.  [finish  ?]  if  x  =  junction  then  go  to  d. 
f.3.  [ratio  test]  the  same  as  l.d.l  -  l.d. 4. 
f.4.  go  to  f.2 

3.  Travel  along  the  loop 

a.  [i  and  j  in  the  same  subtree  ?]  if  R1  »  R2  then  go  to  9* 

b.  [i  and  j  in  the  same  1-tree  ?] 


if  1(R1)  »  1(R2)  f  0  then  go  to  f. 


1,  continued. 

c.  [evaluate  y(x)  along  Rl's  loop[  ] 

c. l  [initialize]  x  :=  R1 

w  :=  R1 
T  :=  T1 

C.2.  [self-loop  ?]  if  l(w)  =  0  then  y(w)  :=  T 

go  to  d. 

C.3.  [loop  factor]  if  a(w)  >  0 

then  T  :=  T/[l  -  factor  (  -  l(w))] 
y(w)  =  y(w)  +  T 
T  :=  -T  •  (g(w)) 
go  to  d. 

otherwise  T  :=  T/ |[l-Factor(-l  (w) )]xg[-a(w)][ 
y(w)  =  y(w)  +  T 
T  :=  -T 
go  to  d. 

C.4.  [finish  ?]  if  p(x)  =  w  then  go  to  e. 

C.5.  x:  a  p(x) 

C.6.  [evaluate  y(x)]  the  same  as  l.c. 

d.  [ratio  test] 

d. l.  the  same  as  l.d.l  -  l.d.4. 

d. 2.  go  to  C.4. 

e.  [evaluate  y(x)  along  R2's  loop  node  and  do  ratio  test] 

e. l.  [initialize]  x  :*  R2 


e.2.  the  same  as  c.2  -  c.6  ;  d.l,  d.2,  except 
changing  c.5  go  to  e  to  c.5  terminate  , 


and  changing  L  =  1  to  L  =  2 
^  ^  u  u 

f.  the  node  i,  j  in  the  same  1-tree  but  have  different  root  node 
[start  from  R1  travel  the  loop,  but  only  y(x)  is  evaluated] 
the  same  as  c.l.  -  c.6.  except  changing  "go  to  d"  to 

"go  to  C.4" 

g.  [combine  two  travel  as  one] 
g.l.  [initialize]  x  R2 

w  :=  R2 
T  :=  Tl  +  T2 
g.2.  the  same  as  e.2. 


In  order  to  obtain  the  proper  dual  variables  and  level  indices,  a  sub¬ 
tree  needs  to  be  traversed..  This  traversal  is  effected  by  making  use  of  the 
successor  and  brother  labels.  To  traverse  the  entire  subtree,  the  following 
rule  is  employed: 

(a)  If  the  node  has  a  successor,  visit  it;  otherwise 

(b)  If  the  node  has  a  brother,  visit  the  brother,  otherwise 

(c)  If  the  predecessor  of  the  node  is  the  root  of  the 
subtree,  terminate;  otherwise,  visit  it. 

The  above  is  the  essential  traversal  which  is  used.  The  particular 
pivot  operation  and  the  dual  variable  update  dependent  on  the  particular 
type  of  pivot  to  be  performed.  If  the  entering  arc  is  in  the  same  1-tree 
of  the  forest  then  the  operation  is  a  little  more  complex  than  the  case 
in  which  the  two  incident  nodes  of  the  entering  arc  are  in  different  1-trees. 

At  the  termination  of  the  ratio  test  algorithm,  the  maximum  allowable 
change  A  is  given,  and  if  u  =  0,  then  the  entering  arc  is  also  the  leaving 
arc.  In  this  case  only  a  flow  change  is  required.  Otherwise,  the  leaving 
arc  is  u,  which  is  on  the  path  L^.  The  path  from  the  tail  node  of  the  entering 
arc  to  its  root  node  is  considered  to  be  path  1  and  the  other  path  as  path  2. 

The  various  possibilities  in  terms  of  the  type  of  change  in  the  basis 
occurring  due  to  the  specifics  of  the  entering  arc  and  the  leaving  arc  have 
been  sumnarized  in  the  introduction.  The  whole  operation  is  separated  into 
three  parts  for  exposition:  Flow  Update,  Forest  Update  and  the  Dual  Variable 
Update. 

The  flow  update  is  rather  straightforward  and  is  summarized  in 


Algorithm  2. 


Algorithm  2.  Flow  Change 

1.  [have  junction  node?]  if  junction  =  0  then  go  to  5 

2.  [update  flow  along  path  from  i  to  junction] 

a.  [initialize]  x:  =  i 

b.  [finish?]  if  x  =  junction  then  go  to  3 

c.  [update  flow]  f(x):  =  f(x)  -  5  •  A  •  y(x) 

d.  [next  node]  x:  =  p(x)  ,  go  to  b 

3.  [update  flow  from  j  to  its  root] 

a.  [initial ize]  x:  =  j 

b.  [finish?]  if  x  =  R2  ,  then  go  to  4 

c.  [update  flow]  f(x):  =  f(x)  -  o  •  A  •  y(x) 

d.  [next  node]  x:  =  p(x) 

go  to  b. 

« 

4.  [update  flow  along  the  loop] 

a.  [initialize]  x:  =  R2  ,  w:  =  R2 

b.  [update  flow]  f(x):  =  f(x)  -  o  •  A  •  y(x) 

c.  [finish?]  if  p(x)  f  w  then  x:  =  p(x) 

go  to  b. 

d.  terminate 

5.  [update  flow  from  i  to  its  root] 

a.  [initialize]  x:  »  i 

b.  [finish  ?]  if  x  »  R1  then  go  to  6 

c.  [update  flow]  f(x):  »  f(x)  -  5  •  A  •  y(x) 

d.  [next  node]  x:  »  p(x) 

go  to  b. 

6.  [How  many  1  -  tree?]  if  1(R1)  »  1(R2)  f  0  or  R1  »  R2 

then  go  to  3. 


Algorithm  2,  continued. 

7.  [update  flow  along  the  loop  of  the  1 

a.  [initialize]  x:  *  R1  ,  w:  »  R1 

b. c.  the  same  as  4.b  and  4.c 
d.  go  to  3 


-  tree  which  contain  node  i] 


In  order  to  update  the  forest.  Algorithm  3  is  effected.  This 
procedure  makes  use  of  the  information  obtained  during  the  ratio  test 
which  determines  the  kind  of  update  required.  That  is,  the  leaving 
arc  is  associated  with  node  u  and  indicates  the  path  which 

contains  node  u*  The  entering  arc  k  is  from  node  i  to  node  j. 
Algorithm  3.  Update  Predecessor  orientation 


1. 

[Initialize]  if  L,  =  1  then  x: 

.mI 

*  i  ,  w: 

=  j 

otherwise 

*  j  ,  w. 

®  i 

2. 

[store  some  old  data  related  to  x] 

Pg*.  ’  P(x) 

T:  »  -a(x) 

f^:  ’  T{x) 

3. 

[delete  x  from  its  old  family] 

a. 

[is  X  the  successor  of  p^^?] 

if  X  »  s(Pq)  then 

s(p„): 

*  b( 

go  to 

4 

otherwi se 

z:  •  s(pg) 

b. 

[Is  X  thfr  brother  of  z?] 

if  X  »  b(z)  then 

b(z):  ■ 

b(x) 

go  to 

4 

c. 

[check  next  one]  z:  *  b(z)  , 

go  to  b 

• 

4. 

[set 

w  be  the  predecessor  of  x] 

p(x):  ■ 

w 

Algorithm  3,  continued. 


5.  [set  X  be  the  successor  of  y] 

b(x) : 

=  s(w) 

s(w): 

»  X 

f(x): 

»  A 

a(x) : 

6.  [finish?]  if  x  =  u  then  terminate 

7.  [next]  w:  =  X 


go  to  2. 


Note  that  if  the  leaving  arc  is  a  loop  arc,  the  pivot  will  break 
the  old  loop.  For  such-  a  pivot.  Algorithm  3  must  be  preceded  by 
making  the  loop  node  a  successor  of  its  predecessor.  The  update  of 
the  forest  is  summarized  in  Algorithm  4. 

Algorithm  4.  Forest  Update 

f 

1.  [set  flag  for  new  loop]  New:  »  0 

2.  [how- many  L  -  tree?]  if  1(R1)  »  1(R2)  j*  0  or  R1  »  R2 

then  go  to  4- 

3.  [two  1  -  tree  case,  no  new  loop  formed] 

a.  [old  loop  broken?]  if  Uu)  >0  then  go  to  c 

b.  [make  the  loop  node  be  the  successor  of  its  predecessor] 

b.l.  [initialize]  x:  •  u 

w:  *  p(x) 

b.2.  [make  x  be  the  w's  successor] 

b(x):  ■  s(w)  , 
s(w):  •  X 


%■  v'^ 


-r.-AN'-'-V-V - 
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Algorithm  4,  continued. 


b.3.  [finish?]  if  w  »  u  then  go  to  c 


otherwise  x:  =  w 


w:  =  p(w) 
go  to  b.2. 


c.  Use  Algorithm  3. 

d.  terminate. 


4.  Only  one  1  -  tree,  a  new  loop  might  be  formed 

a.  [a  new  loop?] 

if  junction  >  0  and  1  (p  )  ^1  (junction) 
then  New  =  1 

and  if  1  ( i )  <  1  ( j )  then  »  1 

otherwise  L  =  2 
u 

b.  [old  loop  broken?]  if  l(p)  >0  then  go  to  e 

c.  [old  loop  broken,  new  loop  must  be  built]  New  »  1 

d.  [make  the  loop  node  be  the  successor] 

the  same  as  3.b. 
a.  Use  Algorithm  3- 
f.  if  New  »  0  then  terminate 

5.  Delete  the  node  on  the  new  loop  frirni  its  predecessor's  family 

a. ,  [initialire]  x:  *  i  ,  w:  »  p(x) 

b.  [delete  x  from  w's  family] 

b.l.  If  X  *  s(w)  then  s(w):  »  b(x) 

go  to  c 
otherwise  w  *  s(w) 

b.2.  if  X  •  b(w)  then  b(w):  *  b(x) 


go  to  c 


Algorithm  4.  continued. 

b.3.  w:  =  b(w)  ,  go  to  b.2. 

c.  [finish?]  if  i  /  p(x)  then  x:  »  p(x) 

w:  =  p(x) 
go  to  b. 

d.  terminate. 

Specifics  of  the  update  of  the  dual  variables  and  the  level 
index  are  summarized  in  Algorithm  5.  Recall  that  the  dual  variables 
are  the  solution  to  the  following  system  of  equations: 

v'B  =  c' 

The  dual  variables  are  updated  since  the  current  basis  and  the 
updated  basis  differ  in  only  one  column.  Only  the  following  sub¬ 
system  need  be  considered: 


Thus,  Vjjd^  +  *  ^ic 


v_d^  +  v_g_  •  c_ 

q  q  p»q  q 


The  subMtrix 


corresponds  to  the  loop  of  1-tree. 
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If  there  no  new  loops  formed,  the  above  submatrix  remains  unchanged;  therefore 
(Vp,...v^)  satisfy  the  last  p-q+1  equation  of  the  above  system,  and  other 
(Vi,...Vp_i)  can  be  easy  calculated  by  the  recursion  equation  traveling 
up  the  subtree.  Hence  the  key  point  is  when  a  new  loop  is  formed,  we 


Algorithm  5,  continued. 

2.  at  least  one  of  node  i  or  j  has  correct  dual  variable  and  level  index 

a.  [initialize]  if  ^  then  w  :=  j  ,  x  :=  j 

b.  [check  successor]  d  :=  x 

X  :=  s(z),  if  X  =  0  then  go  to  e. 

c.  [direction  ?]  if  a(x)  >0,  then  v(x)  ;=  c(a(x))  -  g[a(x)]  •  v( 

otherwise  v(x)  :=  [c(-a(x))  -  v(z)]/g[-a(x)] 

d.  [update  level]  l(x)  =  l(z)  +1  go  to  b. 

e.  [finish  ?]  if  z  =  w  then  go  to  g. 

f.  [check  brother]  x  :=  b{z)  ,  z  =  p(z) 

if  X  0  then  go  to  c. 
otherwise  go  to  e. 

g.  terminate. 

3.  evaluate  the  dual  variable  of  new  loop  node  i,  assume  the  name  of 
new  loop  is  a  . 

a.  [initialize]  x  :=  i 

Vo  :=0 

T  :=  1 
F  :=  1 

b.  if  z(x)  >  0  then  Vq  Vo  +  T  •  c(x) 

T  :=  -T  »  g(x) 

F  :=  -g(x)  .  F 
Otherwise  T  T/g(x) 

Vo  ^0  +  T  •  c(x) 

T  :=  -T 
F  :»  -F/g(x) 

c.  [finish  ?]  if  p(x)  *  i  then  go  to  d. 

otherwise  x  :»  d(x)  .  go  to  b. 


Algorithm  5,  continued 

d.  [store  loop  vector]  Factor  (a)  :=  F 

e.  [store  v(i)]  v(1)  :=  v^/d-F) 

4.  update  dual  variable  and  level  by  traveling  entire  new  1-tree. 

a.  [initialize]  x  ;=  i  go  to  c. 

b.  [set  dual  variable  of  the  loop  node  x] 

if  a(w)  >  0  then  v(x)  :=  [c(a(w))  -  v(w) ]/g[a(w) ] 
otherwise  v(x)  :=  c(a(w))  -  v(w)  •  g(-a(w)) 

c.  l(x)  :=  0  ,  w  :=  x. 

d.  [travel  the  subtree  rooted  at  node  x] 

the  same  as  2.b.  -  2.f. 

g.  [set  level  index  for  node  w]  l(w)  :=  -a 

h.  [finish  ?]  if  p  (w)  f  i  then  x  :=  p(w)  ,  go  to  b. 


otherwise  terminate. 


V.  Computational  Results 


The  data  structures  and  algorithms  presented  in  earlier  sections  have 
been  implemented  in  FORTRAN  to  develop  a  generalized  network  code  at  the 
Center  for  Cybernetic  Studies,  The  University  of  Texas  at  Austin.  The  com¬ 
putational  testing  for  this  code  has  been  carried  out  on  two  sets  of  pro¬ 
blems:  Pure  network  problems  whose  specifications  are  given  in  Table  1,  and 
Generalized  Network  problems  whose  specifications  are  given  in  Table  3. 

These  problems  have  been  generated  using  NETGEN  [  5  ]  and  one  of  its  variants. 
The  pure  network  problems  were  solved  using  a  pure  network  code;  this  code  and 
the  generalized  network  code  were  both  developed  at  the  Center  for  Cybernetic 
Studies.  The  details  of  this  computational  testing  are  given  in  Table  2. 

This  testing  was  carried  out  in  order  to  determine  the  overhead  required  by 
the  generalized  network  data  structures  over  pure  network  data  structures. 
Degradation  in  solution  time  varies  from  a  factor  of  2.2  to  a  factor  of  4.0 
while  the  degradation  in  time  per  pivot  varies  from  a  factor  of  2.6  to  3.7. 

The  results  indicate  that  generalized  network  algorithms  are  on  the  order  of 
2.5  to  3.5  times  slower  than  pure  network  algorithms. 

The  generalized  network  problems  solved  ranged  from  200  nodes  to  1000 
nodes  and  1500  arcs  to  7000  arcs.  The  timings  obtained  for  the  solution  of 
these  problems  are  given  in  Table  4.  The  results  indicate  that  the  data 
structures  presented  are  well -suited  to  carry  out  primal  simplex  operations 
for  generalized  network  problems. 
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for  efficient  solution  of  generalized  network  problans  necessarily  depends  on 
data  structures  chosen  to  represent  the  basis.  Ihis  paper  presents  the  desi 
and  detailed  algorithmic  specification  of  the  primal  simplex  algorithm  for 
sudi  problans.  CcoDutational  testing  to  determine  the  overhead  required  by 
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